Entanglement Spectra of the quantum hard-square model: 
Holographic minimal models 
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We study the entanglement properties of a quantum lattice-gas model for which we can find 
the exact ground state (of the Rokhsar-Kivelson type). The ground state can be expressed as a 
superposition of states, each of which is characterized by a particle configuration with nearest- 
neighbor exclusion. We show that the reduced density matrix of the model on a ladder is intimately 
related to the transfer matrix of the classical hard-square model. The entanglement spectra of 
the model on square and triangular ladders are critical when parameters are chosen so that the 
corresponding classical hard-square models are critical. A detailed analysis reveals that the critical 
theories for the entanglement Hamiltonians are c < 1 minimal conformal field theories. We further 
show that the entanglement Hamiltonian for the triangular ladder is integrable despite the fact that 
the original quantum lattice-gas model is non-integrable. 
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I. INTRODUCTION 

Recently, there has been a growing interest in under- 
standing the characteristics of quantum entanglement in 
many-body systems [l|,l2|. Entanglement (von Neumann) 
entropy and the family of Renyi entropies have been 
demonstrated to be useful in quantifying bipartite en- 
tanglement in a variety of situations. More recently, the 
concept of entanglement spectrum (ES) introduced by Li 
and Haldane [|| has attracted considerable interest. The 
ES of a bipartite system consisting of A and B is obtained 
from the spectrum of the reduced density matrix pa for a 
subsystem A. Li and Haldane demonstrated that the en- 
tanglement spectra for fractional quantum Hall states re- 
flect the gapless edge excitations at the fictitious bound- 
aries created by tracing out the subsystem B. Since then, 
entanglement spectra have been studied in quantum Hall 
systems [341], topological insulators (lp| - [l3l. and quan- 
tum spin models in one [1 3l — ll 8j j and two [19|-|22| dimen- 
sions. 

In Ref. [lfj, Poilblanc studied the ES of gapped two- 
leg spin- 1/2 Heisenberg ladders and found that the ES is 
remarkably close to the energy spectrum of a single spin- 
1/2 Heisenberg chain, which is one of the subsystems of 
the original two-leg ladder. This observation was later 
explained analytically using first-order perturbation the- 
ory in the limit of strong rung coupling [l8|, [23| . It is 
natural to expect that a similar property also holds for 
gapped (but non-topological) spin systems in two dimen- 
sions since the entanglement in those systems is short- 
ranged. From this perspective, the two-dimensional 
(2D) valence-bond-solid (VBS) states [22j |. which are 
exact ground states of the Afflcck-Kennedy-Lieb-Tasaki 
(AKLT) model [U,[2^|, and their generalizations called 
projected entangled pair states (PEPS) [20j were exam- 



ined. In both cases, it was found that the reduced density 
matrix pa can be interpreted as a thermal density matrix 
of a holographic (fictitious) one-dimensional (ID) system. 
In particular, for the VBS states, this holographic ID 
system turns out to be a spin-1/2 Heisenberg chain. 

In this paper, we introduce and study another class of 
quantum many-body states, which falls into the category 
of tensor- network states |26l - l30| . The quantum lattice- 
gas model in which these states are ground states was first 
constructed by Lesanovsky [3l[ . Similar to the Rokhsar- 
Kivelson (RK) point in the quantum dimer model [32j . 
the ground states are expressed as a weighted superposi- 
tion of states labeled by classical configurations of parti- 
cles. All the allowed configurations must satisfy the con- 
straint that there is no more than one particle at any pair 
of neighboring lattice sites. The weight of each configu- 
ration depends on the parameter z, which corresponds to 
the activity (fugacity) in the classical lattice-gas model. 
The parent Hamiltonian of this state in one dimension 
was considered in Ref. [3l[ . It was shown that the ground 
state of the ID Rydberg lattice gas can be well approx- 
imated by this state. The model on the square lattice, 
referred to as the quantum hard-square model, was also 
examined for the validity of this ansatz 33]. Interest- 
ingly, it was found that the normalization of the ground 
state in this case is given by the partition function of the 
classical hard-square model 34] . 

While the parent Hamiltonian can be easily con- 
structed for any graph in any dimension, we mainly fo- 
cus on the model on two-leg ladders in this paper. In 
this case, the ES can be obtained from the eigenvalues 
of the transfer matrix in the classical lattice gas model 
in two dimension. From this, it immediately follows that 
the ES is gapless when the parameter z is chosen so that 
the corresponding classical model is critical [43[ . We nu- 
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merically study the ES of our models on both square 
and triangular ladders at the critical point and find that 
the low-energy part of the entanglement Hamiltonian for 
the case of square (triangular) ladder is well described 
by the conformal field theory (CFT) with central charge 
c = 1/2 (c = 4/5). To our knowledge, this is the first 
example of quantum many-body systems whose entan- 
glement Hamiltonians are precisely described by c < 1 
minimal CFTs. This conclusion is further supported by 
the nested entanglement entropy (NEE), the von Neu- 
mann entropy of the ground state of the entanglement 
Hamiltonian. For the model on the triangular ladder, we 
show that the entanglement Hamiltonian is integrable in 
the same sense as in the hard-hexagon model solved by 
Baxter (jy). 

The organization of the rest of the paper is as follows. 
In Sec. II, we define the model considered and construct 
the exact ground state. We then show that the normal- 
ization constant of this state is closely related to the par- 
tition function of the classical lattice-gas model in the 
same dimension. In Sec. Ill, we consider our model on 
two-leg square and triangular ladders. We show how to 
obtain the spectrum of the reduced density matrix of the 
system from the spectrum of the two-row transfer matrix 
in the corresponding 2D classical model. Then numerical 
results of the entanglement entropy (EE) and correlation 
length are shown. The ES and the NEE at the critical 
point are studied in detail in the same section. We con- 
clude with a summary in Sec. IV. In Appendix A, we 
give a detailed proof that the constructed ground state 
has zero-energy. In Appendix B, we prove the integra- 
bility of the reduced density matrix for the model on the 
triangular ladder. 



II. THE MODEL 

A. Hamiltonian 

The model we introduce is a slight generalization of the 
model constructed by Lesanovsky [3l| . It can be defined 
on any lattice in any dimension. The model describes 
interacting hard-core bosons on a lattice. The Hilbert 
space at each site is spanned by \rii), where m = 1 (0) 
indicates that the site i is occupied (empty). Note that 
there is at most one particle at each site. We further 
assume that the entire Hilbert space is restricted so that 
there is no more than one boson on any pair of nearest- 
neighboring sites. With the identification | f) f-> |1) and 
| I) <-> |0), the operator that creates a hard-core boson 
at site i is expressed as a +V {l} , where a ± = (a x ± ia v )/2 
with a a (a = x, y, z) being the Pauli matrices. The pro- 
jection operator Vus requires all sites adjacent to site i 
to be empty. More explicitly, 

V {l) =\[{l-n 3 ), (1) 



where rij := rf&j = (f| + l)/2 and Gi is a set of sites 
adjacent to site i. It is easy to see that (P{i)) 2 = "Pus- 
Note that commutes with any operator at site i. 

Let A be a lattice (graph). The Hamiltonian of the 
model on A is given by 

H = Y,h\{z)hi{z), h i (z) = [ar-^(l-n i )]r {i} . (2) 

We assume that the parameter z is real and nonnegative. 
We will see later that z exactly corresponds to the activ- 
ity of the classical hard-core lattice gas model. To see the 
physical meaning of the Hamiltonian, it is illuminating 
to rewrite H as 

H = - VI^>+ + ar)V (l) + - z)n t + z}V {i) . (3) 

The first term corresponds to the creation and anni- 
hilation of hard-core bosons, while the second term 
includes the chemical potential term and interactions 
among bosons that are more than two sites apart. For a 
ID periodic chain of length L, the Hamiltonian reads 

L 

H = Y, P[-Vx>i + (1 - + zm-in i+1 + z]V, (4) 

where V denotes the projection operator onto the physi- 
cal subspace in which any pair of neighboring sites can- 
not be occupied simultaneously. One can interpret this 
Hamiltonian as the quantum Ising chain in a trans- 
verse and longitudinal field with nearest-neighbor ex- 
clusion and further neighbor interactions. In higher di- 
mensional lattices, the Hamiltonian involves complicated 
and presumably unphysical multi-body interactions, e.g. 
HiTijUk, and it cannot be written as the quantum Ising 
model except in special lattices such as the 3-12 (Fisher) 
lattice. 

B. Exact ground state 

The Hamiltonian Eq. ([3]) in one dimension and that 
in 2D square lattice were studied in Refs. [3l|, H3[. It 
was shown that the exact ground state of H can be ob- 
tained analytically. In addition, this exact ground state 
resembles the ground state of the Rydberg lattice gas 
which is experimentally realizable in cold atom systems. 
In the ID model, explicit expressions for the energy gap 
and a couple of excited states were also obtained analyt- 
ically m. 

Here we consider the model on a generic lattice (graph) 
and find the exact ground state. Since each local Hamil- 
tonian h\(z)hi(z) is positive semi-definite, the Hamilto- 
nian H is also positive semi-definite and thus the energy 
eigenvalues are nonnegative. Therefore, a state with zero 
energy is a ground state and is annihilated by all hi(z). 
One can find that the zero-energy state takes the form 

I*) = -^== II exp(V^^ + ^<,))l 4), (5) 
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where S(z) is the normalization constant and | JJ.) 
denotes the all-down state, i.e., | || ... \), corre- 
sponding to the vacuum state. Note that the or- 
der of the product in Eq. J5]) is arbitrary because 
[exp(y/za^V^)),exp(y/zaj'V^))] = for all For a 
detailed proof of H\z) — 0, see Appendix [A] We can 
further show that the zero-energy state \z) is the unique 
ground state of H . This can be seen by noting that all 
of the off-diagonal elements of H are nonpositive and 
satisfy the connectivity condition, and thus the Perron- 
Frobenius theorem is applicable. 

Similar to the AKLT model [37[ and the RK model [H| 
as well as their generalizations called PEPS (38|, the 
ground state \z) can be expressed as a weighted super- 
position of allowed states. To see this, it is instru ctive to 
write the unnormalized ground state | \I/ (z)) = y/E,(z)\z) 
in terms of classical configurations: 

\*(z)) = J2z nc/2 \C), (6) 
ces 

where C label classical configurations of particles on A 
and S is the set of configurations with nearest-neighbor 
exclusion. \C) correspond to the basis states in the quan- 
tum model and they are orthonormal, i.e., (C\C) — Sc,c- 
In the expression Eq. ©, tiq counts the number of bosons 
in the state \C). It is now clear that the normalization 
constant S(z) is the partition function of the classical 
hard-core lattice gas model: 

hooh*(*)I*(*0> = (7) 

ces 

and the parameter z can be interpreted as an activity 
(fugacity) of the lattice gas. This classical model in two 
dimensions has been studied extensively in the literature 
of statistical mechanics. For instance, the model on the 
square lattice is known as the hard-square model. One 
can easily imagine that the system undergoes a phase 
transition from a liquid state at z <C 1 to a solid state at 
z>l, where the difference of the sublattice occupations 
is nonzero. In fact, previous studies have shown that the 
model exhibits an order-disorder phase transition at z = 
z c ~ 3.796 [39i - l4TI | . This second-order phase transition 
belongs to the same universality class as the 2D Ising 
model. The model on the triangular lattice is called the 
hard-hexagon model, which is integrable and was solved 
by Baxter [341 ] . It also exhibits an order-disorder phase 
transition at z = z c = (11 + 5V5)/2 = 11.09017..., but 
this transition belongs to the universality class of the 
three-state Potts model. 

Let us return to the original quantum model. A par- 
ticular feature of the present construction is that correla- 
tion functions in the quantum ground state are the same 
as those of the corresponding classical model |36j . This 
provides useful information on the ground state phase di- 
agram of the quantum model. In fact, if the parameter 
z is tuned so that the corresponding classical model is 
critical, then the ground state has algebraically decaying 
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FIG. 1: (color online) (a) Two-leg square ladder and (b) 
two-leg triangular (zigzag) ladder. The systems are divided 
into two parts, A and B, as indicated by the dashed lines. 
We impose periodic boundary conditions in the leg direction: 
T L +i = Ti and (Tl+i = Ti. 

correlation functions, which suggests that the quantum 
model is critical as well. This type of quantum critical 
point is called conformal quantum critical point because 
the ground state wave functional is conformally invariant 
in the scaling limit [42j . 

III. ENTANGLEMENT PROPERTIES 

In this section, we consider the exact ground state 
Eq. O on a two-leg ladder as shown in Fig. [1] The 
classical partition function Eq. (J7J is defined on the same 
ladder, which is still ID. This implies that there is no crit- 
ical point at finite z and the correlation functions of local 
operators in the ground state are exponentially decay- 
ing. Nevertheless, we find an intimate relation to CFTs, 
which will be revealed by the analysis of the entanglement 
properties. We show that the ES of this system can be 
inferred from the spectrum of the transfer matrix in the 
corresponding 2D classical system. We also study the 
von Neumann entropy associated with the ground state 
of the entanglement Hamiltonian, and show that the un- 
derlying CFT is a unitary minimal model with central 
charge c < 1. 

A. Reduced density matrix 

Consider the two-leg ladder of length L with periodic 
boundary conditions in the leg direction. We divide the 
system into two subsystems, A and B as shown in Fig. 
[T] Let r = {ti, T2, t£\ be a particle configuration 
on the chain in A, and a = {o~x, <72> •■■■> <?l} be that on 
the chain in B. Here r, = 1 (0) if the site i on A is 
occupied (empty). The same applies to Oj. Let |r) and 
|cr) denote the corresponding basis states in the quantum 
model. Then the (unnormalized) ground state |^(z)) on 
the ladder is written as 

r a 
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(a) d c 

= w(a, b, c, d) 

1 z 1 ' 4 z 1 ' 4 z 1 ' 4 z 1 ' 4 z 1 ' 2 z 1 ' 2 

1 z 1 ' 4 z 1 ' 4 z 1 ' 4 z 1 ' 4 z 1 ' 2 

FIG. 2: (color online) (a) Graphical representation of the 
local Boltzmann weights. Allowed configurations and their 
Boltzmann weights for the square ladder (b) and the triangu- 
lar ladder (c). 

where 

L 

[T(z)] T>a := JJu)(tTi,(Ti + l,Ti + l,Ti), (9) 

t=i 

with w(a,b,c,d) being the Boltzmann weights for each 
face (see Fig. [5]). More explicitly, for the square ladder, 
we have 

L 

[T(z)]„ - J]z^+^'/ 2 (l-a l r 4 )(l- ( T 4 a l+1 )(I-T 4 r. i+1 ), 
i=l 

(10) 

and for the triangular ladder, we have 

L 

[t(z)] t ,, = n^ +T,)/2 (i-^«)(i-w+i) 

i=l 

X (l-Ti7-i + l)(l-Ti<7i + i). (II) 

One can think of T(z) in Eq. © as an A^-dimensional 
matrix, where Nl is the number of allowed configura- 
tions in each chain 35]. In fact, we can identify T(z) as 
the transfer matrix of the 2D classical lattice gas model 
with hard-core exclusion [34j. Note that a similar trans- 
fer matrix formalism was also applied to entanglement 
entropies of the 2D RK wavefunction [44| . 

The reduced density matrix for A describing the en- 
tanglement between the two subsystems is defined by 
Pa '■= Tr b[\z)(z\]. To obtain the spectrum of pa, we fol- 
low the approach used in Ref. [45j. A similar approach 
has been applied to the 2D VBS states [22l l4^|. We first 
write the state \z) as 

\z) = ^L=y.\^)®\°)' (12) 



states in B, i.e., \a), are orthonormal, one can easily trace 
out the degrees of freedom in B and obtain 

PA = 7=rrT^\<P<r){<P<,\- (13) 
-{z> a 

We next introduce the Gram matrix M (overlap matrix) 
whose matrix elements are given by 

AW = ^(wM- (14) 

Since the states \t) (in A) are also orthonormal, one finds 

M=^L[T(z)] T T(z), (15) 

where the superscript T denotes matrix transpose. From 
the argument in [45j, we can show the following prop- 
erties: (i) Tr[M] = 1, (ii) all the eigenvalues of M are 
nonnegative, and (iii) nonzero eigenvalues of M and pa 
are identical. It follows from (iii) that the ES associated 
with pa is identical to the eigenvalue spectrum of the 
matrix M. The entanglement Hamiltonian of our model 
is then defined via M := exp(— H~e). One can regard M 
as a thermal density matrix of an auxiliary ID model, 
which we call the holographic model. 

Let us now discuss the relation between the spectrum 
of M and that of the transfer matrix T(z). For the square 
ladder, we have wfa, a i+ i, r i+ i, n) = w(n, r i+ i, a l+1 , a l ) 
which yields [T(z)] T = T(z). We thus obtain M = 
\T(z)} 2 /S(z). The matrix {T(z)\ 2 can be interpreted as 
a two-row transfer matrix by which the classical hard- 
square model is described (see Fig. [3] (a)). For the 
triangular ladder, the local Boltzmann weights do not 
have the above symmetry. Instead, we can interpret 
the matrix [T(z)] T T(z) as a product of transfer matri- 
ces defined on the three- leg ladder shown in Fig. [3] (b). 
One can view the lattice as a triangular lattice deformed 
into a topologically equivalent square lattice with diago- 
nal bonds. Therefore, the two-dimensinal classical model 
corresponding to our quantum model on the triangular 
ladder is the hard-hexagon model. For both square and 
triangular ladders, the entanglement spectra can be read 
off from the spectra of transfer matrices [47|. From this 
correspondence, it is clear that the ES is critical if the 
activity z is chosen so that the corresponding classical 
model is critical. In Sec. IIII CI we study numerically the 
spectrum of M and confirm that this is indeed the case. 

B. Entanglement entropy (EE) 

We now study the EE which can be obtained from 

S = -Tr[M In M] = -J2 Pa lnPa ' ( 16 ) 



where \tp a ) :— J2 T [T(z)] T _ a \T) . It should be noted that 
the basis states \ip a ) are not orthonormal. Since the basis 



where p a (a — 1, 2, • • • , Nl) are the eigenvalues of M. 
We calculate EE using exact diagonalization for the 
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FIG. 3: Two-row transfer matrix for the square ladder (a), 
and that for the triangular ladder (b). 



square and triangular ladders with periodic boundary 
conditions (see Fig. [IJ. In order to avoid boundary con- 
ditions incompatible with the modulation of the dens- 
est packings, we assume that L = 2 to and L = 3 to 
(to G N) for the square and triangular ladders, respec- 
tively. The dependence of EE on the activity z for 
both models is shown in Figs. 0] (a) and 2] (b). In 
the limit z — > oo, the EEs for square and triangular 
ladders become In 2 and In 3, respectively, irrespective 
of the system size L. This can be understood as fol- 
lows: for large z, the ground state is approximately given 
by a superposition of the ordered states with the max- 
imum density of particles. These ordered states are re- 
lated to each other by translations. Therefore, we have 
PA ~ ^=(|0101...)(0101...| + |1010...)(1010...|) for the 
square ladder and a similar one with period 3 for the 
triangular ladder. As a result, we obtain the observed 
saturation values of EE. In the opposite limit z — > 0, the 
EEs become zero since the ground state is the vacuum 
state (see Eq. ©). In the intermediate region between 
these two limits, the EE shows a non-monotonic depen- 
dence on z and has a peak in both square and triangular 
ladders as shown in Figs. [4] (a) and 3] (b). In both cases, 
the peak position is at about z = z c , i.e., the critical ac- 
tivity in the corresponding classical model, and remains 
almost unchanged with increasing L. 

Let us focus on the entanglement properties of the 
model at z = z c . Figures|U(c) and|H(d) show the scaling 
of the EE S(L) for both the square and triangular lad- 
ders. From these plots, it is clear that the EEs at z = z c 
scale linearly with the system size L (corresponding to 
the length of the boundary between A and B) and thus 
obey the area law: 



S(L) = aL + S 



(17) 



where a and So are the fitting parameters inde- 
pendent of L. For square and triangular lad- 
ders, (a, S ) = (0.2272(3), -0.036(6)) and (a,S ) = 
(0.4001(3), 0.020(5)), respectively. In both cases, S is 
nearly zero, su ggesting that the topological EE intro- 
duced in Refs. [481 |49( is zero in our system. This is 




FIG. 4: (color online) (a) EE (S) of the state Eq. © on 
the square ladder as a function of activity z and (b) the same 
for the triangular ladder. As indicated by the arrows, S in- 
creases linearly with increasing system size L in both cases. 
The dotted horizontal lines indicate, respectively, In 2 and In 3 
which are the EEs at z = oo. (c) Size-dependence of S at the 
critical activity z = z c for the square ladder with different 
sizes L = 4 — 22 and (d) the same for the triangular ladder 
with L = 3 — 21. The dotted lines are least squares fits to the 
last four data points using Eq. (|17|) . The fitting parameters 
are shown in the figures, (e) Correlation length divided by 
system size £(z)/L versus z for the square ladder and (f) that 
for the triangular ladder. 



consistent with the fact that CFTs describing the entan- 
glement Hamiltonians of our system are non-chiral as we 
will see later. 

As discussed in the previous subsection, the Gram ma- 
trix M can be interpreted as a transfer matrix in the cor- 
responding 2D classical model. Therefore, we expect that 
the anomalous behavior of EE at z — z c is attributed to 
the phase transition in the classical model. To study the 
nature of the phase transitions in both square and tri- 
angular cases, we perform a finite-size scaling analysis of 
the correlation length. The correlation length is defined 
in terms of the entanglement gap as 



l 



\n[pW{z)lpW{z)X 



(18) 



where and p^(z) are the largest and the second- 

largest eigenvalues of M at z, respectively. Figures H] 
(e) and 0] (f ) show the correlation length divided by the 
system size L as a function of z for both square and 
triangular cases. Clearly, the curves for different system 
sizes cross at the same point z = z c , which implies that 
the dynamical critical exponent is given by 1. Near the 
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FIG. 5: (color online) Finite-size scaling plot of £(z)/L for 
the square ladder (a) and that for the triangular ladder (b). 
The correlation length exponents of the square and triangular 
ladders are v = 1 and v = 5/6, respectively. 



critical point, we expect that £(z)/L obeys the scaling 
relation 



t(z)/L = f((z-z e )L 1 '»), 



(19) 



where v is the correlation length exponent and /(•) is a 
scaling function. Figures [5] (a) and [5] (b) show plots of 
£(z)/L versus the scaling variable (z — z c )L x / v . For the 
square ladder, we find a good data collapse with v = 1, 
which agrees with the correlation length exponent of the 
2D Ising model (see Fig. [5] (a)). On the other hand, for 
the triangular ladder, the exact value v = 5/6 can be ob- 
tained by noting that the two-row transfer matrix shown 
in Fig. [3] (b) is exactly equivalent to that of the hard- 
hexagon model |50| . We obtain an excellent data collapse 
as shown in Fig. E](b). Note that the exponent v = 5/6 
coincides with that of the three-state Potts model 1511] . 



C. Entanglement spectrum (ES) 



quantum Hall states, the underlying CFTs are non-chiral 
in our case, i.e., there are both left- and right- moving 
modes. In general, the excitation energies from which 
towers are generated have the form 



A Q — Ao — —j—\ n L,a 



hR,a) 



(20) 



where v is the velocity, and hjj tCt and hu^ a are (holomor- 
phic and anti-holomorphic) conformal weights. The sum 
hh,a + h,R,a is called scaling dimension. tiL t a ~ hii.a is 
related to the total momentum k. Comparing the en- 
ergy spectrum obtained numerically with the above re- 
lation, we can identify conformal weights of low-lying 
states. Scaling dimensions of several low-lying states are 
indicated in Fig. [5] From the obtained scaling dimen- 
sions, we conclude that low-energy spectrum of He for 
the square ladder model at z = Zg is described by the 
CFT with central charge c = 1/2 [52|, whereas that for 
the triangular ladder model is described by the c = 4/5 
CFT. The former CFT is identical to that of the 2D crit- 
ical Ising model, while the latter describes the universal- 
ity class of the critical three-state Potts model. These 
are all consistent with the analysis of correlation length 
exponent in the previous subsection. 

We note that the exact scaling dimensions for the tri- 
angular ladder model at the critical point can be ob- 
tained analytically by exploiting the integrability of the 
hard-hexagon model [531 ] . Here integrability means the 
existence of a family of commuting transfer matrices. In 
fact, for the model on the triangular ladder, we can show 
that there is a one-parameter family of commuting ma- 
trices including the Gram matrix M in Eq. (JTSJ as a 
special case (see Appendix [B] for the proof). This prop- 
erty holds even away from the critical point. Therefore, 
the entanglement Hamiltonian of this system at arbitrary 
z is integrable in the same sense as in the original hard- 
hexagon model. 



The anomalous behavior of the correlation length at 
z = z c suggests that the entanglement gap at this point 
vanishes linearly with 1/L. To further elucidate the 
gapless nature of the entanglement Hamiltonian He := 
— InM, we calculate the excitation spectrum of He at 
the critical point. The ES {\ a } a =i,...,N L can be obtained 
from the relation X a = — lnp a . Each eigenstate is la- 
beled by the total momentum k due to the translational 
symmetry in the leg direction. Fig. |H] shows the ES of 
both square and triangular ladders at z = z c . In both 
cases, there is the minimum eigenvalue Ao(:= min Q A Q ) 
at k = 0, and the ES is symmetric about k — and k = tt 
(mod 27r). For the square ladder, the gapless modes at 
momenta k = and k — n are clearly visible. On the 
other hand, for the triangular ladder, they are at k = 0, 
2tt/3, and 4tt/3. 

These towers of energy levels can be identified as those 
of CFTs describing the low-energy spectra of He- As op- 
posed to topologically ordered systems such as fractional 



D. Nested entanglement entropy (NEE) 

To provide further evidence that the entanglement 
Hamiltonians at z — z c are described by minimal CFTs 
with central charge c < 1, we next consider NEE which 
was first introduced in Ref. [22[ . The NEE is the en- 
tanglement (von Neumann) entropy associated with the 
ground state of He- One might think that the scaling 
analysis presented in the previous subsection suffices to 
find the central charges. However, the NEE provides a 
simpler and more clear-cut approach. As discussed in 
Sec. IIII Al the Gram matrix M = exp(— He) can be in- 
terpreted as a thermal density matrix of the holographic 
model in one dimension. This implies that He can be 
written as He — P c sH\ 10 \ with i?hoi being the Hamilto- 
nian for the holographic model. Here f3 e s is the ficti- 
tious inverse temperature and is not universal. We are 
interested in i/hoi rather than He- However, there is no 
unique way to disentangle /3 e ff from i/hoi- To overcome 
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FIG. 7: (color online) NEE for the square ladder (a) and that 
for the triangular ladder (b). Dotted lines are a guide to the 
eye and indicate the slope of c/3, where c = 1/2 and c = 4/5 
are used for the square and triangular ladders, respectively. 
The solid circles indicate NEE for L = 24. The other symbols 
are the same as in Fig. [5] 



FIG. 6: (color online) Low-energy spectra of the entangle- 
ment Hamiltonian He at z — z c for the square ladder (top 
panel) and for the triangular ladder (bottom panel). In both 
cases, the system size is L = 18. The ground state energy 
of He is denoted by Ao. The red open circles indicate posi- 
tions of the primary fields of the corresponding CFTs, while 
the green open squares imply the positions of the descendant 
fields. Discrepancies between the numerical results and the 
CFT predictions are due to finite-size effects. 



this, we first note that the ground state of He is the 
same as that of -ffhoi, corresponding to the eigenvector of 
M associated with the largest eigenvalue. We then recall 
that for ID critical systems, one can read off the underly- 
ing CFT solely from the entanglement entropy obtained 
from the ground state of -ffhoi [H| • Therefore, one can 
obtain the central charge directly from the scaling prop- 
erties of the NEE without evaluating f3 c g or the vel ocity 
v in Eq. (|20l) , as demonstrated in the previous work [22| ■ 

Let us now give the definition of the NEE. We divide 
the system on which the holographic model lives into two 
subsystems: a block of I consecutive sites and the rest 
of the chain. The nested reduced density matrix is then 
defined as 



is expected to behave as 



Tr 



t+i, 



(21) 



where \ipo) denotes the normalized ground state of 
and the trace is taken over the degrees of freedom outside 
the block. Using p(£), the NEE is expressed as 

s(£,L) :=-Tt w [/>(*)kp(0]> ( 22 ) 



where the trace is taken over the states in the block. 

Let us now analyze the numerically obtained NEE in 
detail. Since the low-energy spectra of at z — z c show 
good agreement with those of CFT predictions, the NEE 



s(e,L) = -hx[g(£)] + Sl , 
g {l) = - sin — 

7T \ L 



(23) 
(24) 



where s\ is a non- universal constant [54| . Figure [7] shows 
the NEE for both square and triangular ladders as a func- 
tion of ln[<?(^)]. The data for L — 24 are obtained by the 
power method. The slopes of the dotted lines in Fig.[7]are 
c/3, where c = 1/2 for the square ladder and c = 4/5 for 
the triangular ladder. Clearly, the results are in excellent 
agreement with the formula Eq. (|23p and provide further 
evidence for the holographic minimal models, i.e., the en- 
tanglement Hamiltonians associated with the square and 
triangular ladders at described by unitary 

minimal CFTs with c < 1. 



IV. CONCLUSION 

We have presented a detailed analysis of the entan- 
glement entropy and entanglement spectrum (ES) in the 
ground state of the quantum lattice-gas model on two- 
leg ladders. The exact ground state of the model can be 
obtained as a superposition of states, each of which is 
labeled by a classical configuration with nearest neigh- 
bor exclusion. We have shown that the reduced density 
matrix of one of the legs can be written in terms of the 
transfer matrix in the classical lattice-gas model in two 
dimensions. We then numerically studied the entangle- 
ment properties of the models on both square and tri- 
angular ladders. Both models exhibit critical ES when 
the parameter (activity) z is chosen so that the corre- 
sponding classical model is critical. From the finite-size 
scaling analysis, we found that the critical theory describ- 
ing the gapless ES of the square ladder is the CFT with 
c = 1/2, while that of the triangular ladder is the CFT 
with c = 4/5. This was further confirmed by the analysis 
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of the nested entanglement entropy. We also showed that 
the model on the triangular ladder is integrable at arbi- 
trary z in the sense that there is a one-parameter family 
of matrices commuting with the reduced density matrix. 

It would be interesting to calculate other quantities 
by exploiting the close connection between the entan- 
glement Hamiltonian and the two-dimensional classical 
model. It is likely that the virtual-space transfer matrix 
method [HI] can be applied to the calculation of Renyi 
entropies. Finally, it would also be interesting to ex- 
tend our analysis to the case of truly two-dimensional 
lattices. We note, however, that the model on AT- leg lad- 
ders (N > 2) cannot be treated on the same footing since 
the Gram matrix cannot be written by the transfer ma- 
trix alone. Thus new methods need to be developed to 
analyze two-dimensional systems. 
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Appendix A: Zero-energy ground state 

Here we prove that the state \z) defined in Eq. ([5]) is 
the zero-energy ground state of the Hamiltonian. Since 
the Hamiltonian takes the form H = ^2 ie ^h\(z)hi(z), 
it suffices to see that \z) is annihilated by all hi(z). 
For notational convenience, we introduce the operator 
g\{z) :— exp( y /za^'P^) = 1 + y/zafV^). Recalling that 
[9i( z )i9j{ z ')] = Vi, 3 an d (1 - rii)af = 0, we have 

jeA\{t} 

= n Si(*)i^>- ( Ai ) 

jeA\{i} 

We can further simplify the above expression by noting 
that a~ commutes with gj(z) when i and j are not ad- 
jacent, and a~ g\{z) = when i and j are adjacent. It 
then follows from a7\ JJ.) = that hi(z)\^(z)) = for 
any i € A, which proves that \z) is a zero-energy ground 
state of H. For the uniqueness of the ground state, see 
the discussion in Sec. Ill Bl 



Appendix B: Hard-square model with diagonal 
interactions 

To find a one-parameter family of commuting matri- 
ces for the model on the triangular ladder, we consider 
a larger class of models in which the hard-square and 
hard-hexagon models are included as limiting cases. This 
generalized model was first introduced by Baxter to solve 
the hard- hexagon model [34| . 

Besides the activity z, there are two interaction param- 
eters L and M (not to be confused with the system size 
L and the Gram matrix M as used in the main text) in 
the generalized model. The transfer matrix of the system 
corresponding to Eq. © is defined as 

[T(Z,L,M)] T ^ :=[|lil((7„(J, + l,T 1+1 ,T,) (Bl) 

i 

with 

W(a,b,c,d) = z ("+b+c+d)/i e Lac+Mbd t -a+b-c+d 

x (1 - a6) (1 - be) (l-cd)(l- do) ,(B2) 

where a, 6, c, d take the value or 1. Here we again 
consider the case with periodic boundary conditions in 
the horizontal direction. The parameter t cancels out in 
T(z, L, M) and so is arbitrary. L and M represent the 
diagonal interactions as shown in Fig. [5] (a). The Boltz- 
mann weights for the allowed configurations of each face 
are shown in Fig. [5] (b). The model on the square lad- 
der is obtained by setting L = M = 0, while the model 
on the triangular ladder is obtained by taking the limit 
L = and M -> -co. 

Baxter showed that the statistical mechanics model 
with Boltzmann weights Eq. (|B2j) is exactly solvable if 
z,L,M satisfy 

_ (1 - e~ L ){l - e~ M ) 

Z ~ e L+M _ e L _ e M ' ( B3 ) 

He also showed that two transfer matrices T(z, L, M) and 
T(z', L', M') commute if z, L, M and z', L' , M' satisfy the 
above condition and 

z-^ 2 (l - ze L+M ) = (zT 1/2 (l - z'e L ' +M '). (B4) 

On the integrable manifold represented by Eq. (|B3|) . the 
activity z is determined as a function of L and M . How- 
ever, it is unconstrained in the special limit L = and 
M — > — co, corresponding to the hard-hexagon model. 
Note that the pure hard-square model is not solvable 
since Eq. (|B3[) has only the trivial solution z = if 
L = M = 0. 

Let us return to our model on the two-leg ladder. The 
Gram matrix in Eq. (|15[) can be embedded in the follow- 
ing larger class of matrices: 

M(z,L,M) := l L [T(z,L,M)] T T(z,L,M), 

(B5) 
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(a) d c 




FIG. 8: (a) Graphical representation of the local Boltz- 
mann weights, (b) Allowed configurations and their Boltz- 
mann weights, (c) Row-to-row transfer matrix T(z, L, M) 
and its transpose [T(z, L, M)] T = T(z, M, L). 



where ~(z,L,M) = Tr[T(z, L, M)] T T(z, L, M). Since 
the relations (|B3|) and (|B4|) are symmetric with re- 
spect to the interchange of L and M, we have 
[T(z,L,M),T(z,M,L)] = 0. By the same reasoning, 
[T(z,L,M),T(z',M',L')} = if z',L',M' satisfy both 
Eqs. JB3J and (|Bg|). On the other hand, from Eq. (|B2l) . 
it is easy to see that [T(z,L,M)} T = T(z,M,L) (see 
Fig. [5] (c)). It immediately follows that the matrix 
T(z, L, M) commutes with its transpose. Therefore, four 
matrices T(z,L,M), [T(z,L,M)] T , T(z',L',M'), and 
[T(z' , L' , M')] T are mutually commuting if z,L,M and 
z',L',M' satisfy both the condition (|B3I) and the rela- 
tion (|B4I) . This ensures the existence of the family of 
commuting matrices M{z, L, M). The matrix specializes 
to the Gram matrix for the triangular ladder model when 
L = and M -> -oo. 
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